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NUCLEAR MAGNETIC RESONANCE STUDY COLEMANITE! 


ABSTRACT 


electric temperatures below about —2° has been investigated means 
nuclear magnetic resonance (n.m.r.) techniques. The resonances are split 
because the nuclear Zeeman levels are perturbed the interactions between the 
nuclear electric quadrupole and the electric field gradients existing 
the boron sites. The splittings have been examined detail room temperature 
and —40°C. The results have been analyzed and the quadrupole coupling 
constants, the asymmetry parameters, and the orientations the principal axes 
the electric field gradient tensors existing the boron sites room tempera- 
ture and —40° are reported. Selected resonance lines have been examined 
over the temperature range —120° with particular emphasis the 
region about where phase transition occurs. The complex proton signal 
was also studied over the same temperature range. Abrupt broadening this 
signal occurred the phase transition. These studies revealed that the crystal 
may transform from its centrosymmetrical room-temperature (point group 2/m) 
form either metastable monoclinic form with point group triclinic 
form with point group not clear whether two transitions, separated 
only about are involved whether there only one transition with two 
alternative arrangements, differing only slightly activation energy, available 
the structure. The transition transitions are the second-order displacive 
type. Where possible, the results have been interpreted terms the crystal 
structure. 


INTRODUCTION 


The mineral colemanite, has aroused keen interest 
because was the first boron-containing material known ferroelectric. 
the basis their X-ray investigations colemanite room temperature, 
Christ, Clark, and Evans (1958) proposed centrosymmetrical structure which 
consistent with the space group However, Davisson (1956) observed 
that colemanite weakly piezoelectric and pyroelectric room temperature 
and concluded that colemanite lacks center symmetry. the same time 
Davisson found indications that colemanite undergoes phase transition; the 
pyroelectric constant shows exceedingly high and narrow peak —2.5° 
and anomaly exists the dielectric constant this temperature. Goldsmith 
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(1956) found, agreement with the work Davisson, sharp dielectric 
anomaly centered and also discovered that colemanite ferro- 
electric below this temperature. However, the dielectric breakdown patterns 
(Davisson 1956) remain centrosymmetrical over wide range temperatures 
including that the transition. Furthermore, there optical evidence 
phase change and, the time this work was undertaken, there was X-ray 
evidence either. 

The object this investigation was threefold: (1) obtain quantitative 
data the quadrupole coupling tensors the room temperature and 
ferroelectric phases colemanite, (2) resolve the discrepancies regarding 
the symmetry colemanite, and (3) detect, possible, the phase transition 
and study its progress function temperature. 


THEORY 


nucleus spin and magnetic moment placed magnetic field 
has available equally spaced Zeeman levels. The transitions 
between adjacent levels give rise single absorption line the Larmor 
frequency the nucleus also possesses electric quadrupole 
moment and situated electric field with non-vanishing gradient 
then, Pound (1950) has shown, each the Zeeman levels perturbed, 
and the single line split into components. 

have studied experimentally this so-called quadrupolar splitting the 
3/2) resonance colemanite and have analyzed our data following 
the method developed Volkoff (1953) for the case where the electric quadru- 
pole interaction weak compared with the interaction between the nuclear 
magnetic dipoles and For convenience reference and establish nota- 
tion, summary the relevant theory, specialized the case 3/2, 
reproduced below. Each unique nuclear site gives rise three lines: central 
line from the transition +1/2 —1/2, and two and 
from the transitions +3/2 +1/2. According standard perturba- 
tion theory, the frequencies the central line and the centroid the 
satellites, 1/2 show displacement from provided the quadrupole 
interaction weak that permissible neglect all perturbation terms 
higher than the first. The frequency difference between the two satellites 
given, for rotation the crystal about any axis maintained perpendicular 


where the angle rotation. the case somewhat stronger quadrupole 
coupling, the first-order perturbation term may inadequate and the second- 
order perturbation term must also used. The second-order term contributes 
nothing new the splitting the satellites, but predicts shift from 
the centroid the satellites and the central line. Thus, for rotation about 
any axis perpendicular 


n+p cos sin cos sin 


(3) 
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where stands for the centroid the pair satellites. Note 
that equation (1) neglects perturbation terms third and higher odd order 
while equations (2) and (3) neglect perturbation terms fourth and higher 
even order. 

the crystal rotated turn about each three mutually perpendicular 
axes which are fixed with respect the crystal, then three independent 
6x, and (the subscripts denote the rotation axes). These coefficients are 
related the components the quadrupole coupling tensor 
(where the electrostatic potential the nuclear site and 
expressed terms the axes. For the rotation, the coefficients 
equation (1) are given 


(4) 


Corresponding relations for the rotations about the and axes may 
obtained cyclic permutation the subscripts. consequence 
Laplace’s equation 


These identities are useful check the experimental accuracy. 
The coefficients equation (2) for the shift the centroid the pair 
satellites the rotation are related the coefficients equation (1) 


(5) 


(6) 
dy = — bycx 


Also, equation (3) for the shift the central line the rotation, 
the coefficients can written 


(7) 


Again, corresponding relations for the rotations about the and axes may 
obtained cyclic permutation the subscripts. 

combined with the values obtained directly from the data, 
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lead, through equations (4), values for the components the quadrupole 
coupling tensor the co-ordinate system. The relative, but not 
the absolute, signs these components are determined because this experi- 
ment was not possible find out which the two observed satellites 
corresponded which the two transitions +3/2 +1/2. 

This tensor can then diagonalized standard procedures yield the 
quadrupole coupling constant, and the asymmetry parameter, 
The process diagonalization will also give, with respect 
the axes, the direction cosines the principal axes the 
electric field gradient tensor. For second rank tensor the orientation but 
not the sense the principal axes significant only the relative signs 
the direction cosines are determined. 

far have considered only the interactions the nuclei with and 
with VE. later section, dealing with the proton resonance spectrum, 
shall concerned with magnetic interaction between nuclei. 
From classical point view this may regarded the following manner. 
Each nucleus finds itself not only the externally applied magnetic field Ho, 
but also small local field Hice produced neighboring nuclear dipoles. 
The direction varies from one nuclear site another, depending 
the positions and relative orientations the other nuclear dipoles the 
lattice. Since the field produced magnetic dipole decreases 
nuclear site determined mainly the nearest neighbors. follows that 
the total magnetic field will vary few gauss from one nuclear site 
another. consequence, the resonance condition not sharp and the 
signal will broadened. the interacting nuclei are identical, spin-exchange 
processes will limit the lifetime given spin state. The resultant uncertainty 
the energy the state will cause further broadening the resonance 
signal. 

cases where the nuclei each interact strongly with several neighbors, 
the calculation the shape the resonance signal too complex for 
exact solution except for arrangements high symmetry. Fortunately, Van 
Vleck (1948) has developed quantitative relationship between the nuclear 
configuration and the so-called second moment the line shape. This quantity, 
the average value the square the frequency deviation from the center 
the resonance, defined 


where g(v) normalized function taken represent the shape the 
absorption line plotted versus frequency. The result VanVleck’s theorem 
shows that the second moment can expressed 


where the nuclear magneton; are the nuclear g-factor and spin for 
the resonant nuclei, respectively; are the nuclear g-factors and spins, 
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respectively, for other nuclides the sample; the distance between 
nuclei and the angle between and the vector from nucleus 
and the number nuclei within the volume which the broadening 
interactions are considered confined calculating 

Equation (9) was developed for rigid array nuclei. the event that 
the atoms molecules containing the interacting nuclei may have vibrational, 
rotational, translational motion, the resonance line shape and its second 
moment may profoundly altered. the orientation rate much greater 
than the resonance line-width, the function cos? must averaged 
over the motion resulting decreased line-width and second moment. 

One can now predict the expected behavior the width resonance 
line the temperature varied. very low temperatures the signal from 
the solid material will very broad (of the order several tens kilocycles 
per second) and will remain fairly constant width the temperature 
increased until transition from effectively rigid system rapidly 
reorienting system the signal abruptly narrowed. 


THE COLEMANITE CRYSTAL 
Colemanite, mineral found California, U.S.A., 
usually conjunction with other boron-containing minerals. The crystals are 
generally small, colorless, transparent, prismatic, and elongated along 


The crystallographic data obtained Christ (1953) are summarized below: 


monoclinic; space group: 


cell contents: 

Christ, Clark, and Evans (1958), using X-ray techniques, have established 
the structure colemanite. contains infinite chains running 
parallel the axis, linked together laterally ionic bonds through 
ions form sheets parallel (010). The sheets, turn, are linked together 
through system hydrogen bonds involving the hydroxyl groups the 
chains and the water molecules. The chain element has the composition 
and consists BO; triangle, tetrahedron, and 
tetrahedron linked corners form ring. The structure 
shown projection (001) Fig. 

many ferroelectric materials, hydrogen atoms play dominant role 
the ferroelectric behavior; accordingly, shall pay particular attention 
the positions the hydrogen atoms the colemanite structure proposed 
Christ al. The calcium, boron, and oxygen atoms show prominently 
the electron density projections. However, was necessary for Christ al. 
infer the positions the hydrogen atoms from the configuration these 
heavier atoms. shown Fig. only one oxygen atom the asym- 
metric unit not bonded boron atom; this atom must therefore 
identified part the water molecule. Neighboring water molecules are 
separated 2.73 and, following the arguments Christ al., are joined 
hydrogen bond. Since there center symmetry lying midway between 
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projection the structure colemanite (001). The dashed lines indicate the 
postulated hydrogen bonds. The smallest circles indicate the probable positions the hydrogen 
atoms, with whole circle designating site occupied full time hydrogen atom, and 
half-circle designating site occupied half time hydrogen atom. Information taken from 


Christ al. (1958). 


the two water molecules, must assumed that the hydrogen atom the 
bond contributed half the time one molecule and half the time the 
other. The water molecule lies distances 2.76 and 2.88 from Og, 
respectively. This suggests that normal hydrogen bond links the water 
molecule and and that the water molecule, when not contributing 
hydrogen atom the bond, directing hydrogen bond 
toward The remaining three hydrogen atoms must associated 
hydroxyl groups with and since the remaining four oxygens are 
linked two boron atoms each. The hydroxyl groups and appear 
linked hydrogen bonds Os, since the distances are 
and 2.73 respectively. The hydroxyl group lies from 
the oxygen atom and 2.90 from the oxygen atom the hydroxyl group 
This suggests that directs hydrogen bond toward these two atoms, 
each about half the time. this way, Christ have assigned all the 
hydrogen atoms full-time hydrogen bonds average length 2.74 and 


ORIGIN 
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half-time hydrogen bonds average length 2.87 The existence 
bonds implies disordered state, which, course, may either 
dynamic static. The X-ray data cannot distinguish. 

Christ al. believe that the key the ferroelectric behavior colemanite 
lies the hydrogen bonds. They suggest that the disorder the hydrogen 
bonds dynamic and that, the transition temperature, the hydrogen 
bonds between the water molecules become permanently polarized 
ordered arrangement give rise spontaneous polarization. The 
may also become permanently bonded, below the transition point, 
either 


APPARATUS AND EXPERIMENTAL PROCEDURE 


The nuclear magnetic resonance spectrometer, conventional design, 
consisted oscillating detector the type described Volkoff, Petch, 
and Smellie (1952), narrow-band audio amplifier tuned the magnetic 
field modulation frequency Schuster-type phase-sensitive 
detector (1951), and recording milliammeter. The externally applied mag- 
netic field, produced 12-in. Varian electromagnet, was monitored 
means proton resonance signal displayed oscilloscope. General 
Radio heterodyne frequency meter Type 620-A, with the aid Hallicrafter 
Model SX-62A communications receiver, was used make frequency measure- 
ments with accuracy about 0.5 kc/sec. 

The original crystal mount and dial assembly, which were used for the 
measurements obtained room temperature, permitted the axis, about which 
the crystal was rotated, made perpendicular the magnetic field 
within about minutes arc and the rotation angle read within 
0.3°. For work temperatures other than room temperature, improved 
crystal mount and dial assembly were constructed. The use larger dial 
and circular vernier made possible read the rotation angle 0.1° and 
the arc system from X-ray goniometer incorporated into the crystal mount 
made possible orient the rotation axis that was perpendicular 
the magnetic field within minutes arc. this apparatus, the crystal 
mount, crystal, and coil assembly were all located inside insulated, virtually 
airtight cell which was provided with inlet and outlet. The temperature 
the crystal was controlled blast air which was circulated through 
system heat exchangers series with the cell. was possible maintain 
the temperature the crystal, measured thermocouple, within 
0.5° any desired temperature the range —120° +50° 

number crystals and cleavage fragments colemanite from Inyo and 
San Bernardino counties, California, U.S.A., were used our preliminary 
n.m.r. studies. The individual lines the n.m.r. spectra most these 
samples were found very broad. Ordinarily this would not have been 
too serious, but the case colemanite the spectrum complex with 
large number lines maximum can resolved low temperatures) 
and the separations between many the lines are unusually small that 


J 
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the excessive breadth the signals made impossible resolve them over 
large ranges crystal orientation. Also, some lines were apparently broad- 
ened certain crystal orientations that was not possible detect them. 
Presumably this excessive broadening was due impurities internal stresses 
the crystals. 

For all the work reported here, transparent cleavage fragment cole- 
manite from San Bernardino county, which gave much narrower lines, was 
used. This sample was roughly L-shaped and approximately 
size. possessed two perfect (010) cleavages, one perfect and 
one striated (110) face, and distinct, but imperfect, cleavage parallel 
(001). The intersection, along 2-cm edge, perfect (010) cleavage plane 
and the perfect (110) face defined the axis. The crystal orientation was 
deduced from morphological features and confirmed X-ray diffraction. The 
and axes were chosen parallel the and directions, re- 
spectively, because the external features the crystal made easy recognize 
these directions. The rotation angles and were measured between 
and the and axes respectively. 

For given rotation, the crystal was first mounted accurately possible 
the basis its external features with the desired axis parallel the rotation 
axis and then very accurately oriented observing the the 
spectra. For each the three rotations, the n.m.r. spectra were 
recorded 15° intervals from 360°. Where these data were insufficient 
establish the rotation patterns unambiguously, additional spectra were 
taken intermediate orientations. few orientations certain lines were 
close together that they could not resolved. such cases, the frequencies 
were interpolated from additional readings taken orientations where the 
lines could resolved. 

EXPERIMENTAL RESULTS 
The Spectrum Colemanite Room Temperature 

(a) Qualitative the Spectrum 

three lines, central component and two satellites and would 
expected all the boron sites crystal were equivalent. Actually the 
spectrum colemanite considerably more complex, consisting maximum 
resolved lines shown Fig. 2A. The frequencies the lines 
depend strongly upon the angular orientation the crystal except for that 
the very intense line denoted which remains constant the calculated 
unperturbed value. Whenever the crystal orientation such that the 
crystallographic axis either parallel perpendicular the magnetic field 
direction, number the lines coincide produce spectrum, 
Fig. 2B, seven lines addition the strong central line, vo. The line widths, 
defined the two maxima the derivative each absorption line, were 
the order kc/sec but they did vary some extent depending upon the 
orientation the crystal. The relative intensities did not vary appreciably 
with changes the crystal orientation. 
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K K K K, 


Fic. Selected recordings the derivatives the absorption lines colemanite 
room temperature. The nomenclature the lines explained the text. The frequency scale 
approximately 107 per chart division. 

(A) This recording shows the maximum number resolved lines. was obtained with 
crystal orientation such that the axis made angle 135° with and the and axes 
angles 10° and 80°, respectively, with the normal the plane containing and the 

axis. 

(B) example the reduced spectrum obtained the crystal orientation 150°. 


The frequency variations the lines are plotted Figs. and for 
the three rotations about and Note that for the complete rotation 
the spectrum reduced since the axis maintained perpendicular 
Each pair rotations has one orientation common (e.g., and 

90° both correspond with the axis parallel Ho) lines arising 


K L L K 

[—— —- } 
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Mc/ sec 


Fic. dependence angular position the resonance lines colemanite 
room temperature for the rotation. 


from given nuclear site can identified the various rotations. detailed 
inspection Figs. and reveals that the lines can classified the 
basis their angular dependence. The eight lines denoted 
and are never separated from the frequency the line 
denoted more than about 150 kc/sec and they occur pairs with the 
two members each pair symmetrically situated opposite sides the 


524 
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Mc /sec 
ro) 


11-800 


11-600 


11-400 


ie) 15 30 45 60 75 90 105 120 135 150 165 180 


Fic. dependence angular position the resonance lines colemanite 
room temperature for the rotation. 


line vo, i.e., coincides with vo. They are obviously satellite lines arising from 
four non-identical sets sites which the quadrupole interaction weak. The 
frequencies the central lines arising from these four sites should almost 
exactly equal the unperturbed frequency and independent the crystal 
orientation because the energy difference between nuclear +1/2 and 
—1/2 states virtually unchanged when the quadrupole interaction 
weak. The line much more intense than any other line and, the crystal 
orientation changed, its frequency remains constant the calculated 
unperturbed value. Hence reasonable consider consist the 
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126 


122 


Mc/sec 


11.4 


Fic. The dependence angular position the resonance lines colemanite 
room temperature for the rotation. 


four separate, but unresolved, central signals. The four lines denoted 


M M M M 
and also occur pairs with the two members each 


pair nearly symmetrically situated with respect vo. They must satellite 
signals arising two non-identical sets sites. The quadrupole interactions 
these sites must strong because shifted from vo. Since the quadrupole 
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frequencies the central lines arising from the +1/2 —1/2 
these sites will also shifted from vo. These central lines are clearly the 
remaining two lines denoted and which always remain close 
frequency but not have symmetrically situated (with respect vo) 
counterparts. Thus, after classification, the complex pattern lines reveals 
itself superposition six patterns. 

summary, the lines occurring the n.m.r. spectrum cole- 
manite room temperature can explained one assumes the existence 
six non-identical sets boron sites per unit cell. The electric field gradients, 
which are large two sites and small the other four, must differ, some 
respect, from one set sites another. 

(b) The Point Group Symmetry Colemanite Room Temperature 

Our next aim establish the point group symmetry colemanite 
room temperature. the basis chemical data and the unit cell dimensions, 
colemanite contains boron atoms per unit cell. Remembering that two 
sites unit cell are identical from the point view the n.m.r. experiment 
only they are related center symmetry, apparent that cole- 
manite, room temperature, must have centrosymmetrical structure since 
have found that the boron atoms occupy only six non-identical sites. 
Referring Figs. and can seen that the orientation patterns are 
and that the reduced spectrum occurs the following positions: 
and 90° (Y, that is, parallel Ho), 90° and perpen- 
dicular Ho). Furthermore, appears Fig. reduced spectrum results 
for the complete rotation perpendicular Ho). The reduction dis- 
tinguishable sites three whenever the axis either parallel perpendicular 
the magnetic field can explained only the axis parallel either 
twofold (rotation screw) axis normal mirror (or glide) plane 
both. The combinations these three alternatives with center 
all give the same point group symmetry 2/m. Thus the boron 
sites the unit cell colemanite room temperature may considered 
have been derived the operation the symmetry group 2/m the 
three unique boron sites the asymmetric unit. 

Actually, anticipating our low temperature results, the presence boron 
sites the unit cell colemanite, and hence the center symmetry 
room temperature, can established the analysis n.m.r. spectra alone. 
shall return this point later. 

(c) Nomenclature for the Boron Sites and the Resonance Lines 

The nomenclature the boron sites and resonance lines was adopted with 
full regard paid the fact that the established point group colemanite 
room temperature 2/m. The three sets boron sites are designated 
capital letters and the four symmetry-related sites making each 
set are further identified the subscripts 1,2, and Thus and are 
related the twofold symmetry axis and are and Since and 
(or and are related the center symmetry, these sites 
are identical; hence room temperature, these two sites are simply denoted 
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(or The and sites are similarly denoted. The three lines 
arising from the transitions +3/2 and +1/2 —1/2 are 
denoted and respectively, the subscripts and being chosen 
completed attaching superscript indicate the site its origin, thus 
Each line occurring the reduced spectrum should 
denoted simply the The lines the rotation and 


vo, since frequencies are virtually identical with 


(d) Quantitative Analysis the Spectrum 

For complete rotation the crystal, the frequencies the lines are 
predicted theory periodic with period 180° and fact the 
experimental points and each case were found differ 
very little they were averaged. The Fourier coefficients the curves 
obtained for and were determined the 12-point 
analysis Whittaker and Robinson (1948). The magnitudes the co- 
efficients for pairs, such and which are related the twofold 
screw axis (or the glide plane), were averaged because they possess numerically 
identical VE’s but their relative signs were retained because the VE’s differ 
orientations their principal axes. expected theoretically, the coefficients 
all terms involving angles other than 20, and were found neg- 
ligibly small. The other coefficients, labelled according the nomenclature 
outlined previous section, are listed Table Because their greater 
accuracy, only the coefficients were used the determination the 


TABLE 


Averaged experimental values the Fourier coefficients (in kc/sec) the lines 
colemanite room temperature 


Sites Ki,3, Kas L1,3,L24 Mi,3,Mas 
a +22.5 -7.8 —-11.6 —36.5 —-12.3 +48.4 +641.3  -—675.7 +32.4 
b -6.0 +30.0 —26.7 +63.5 —83.9 +22.6 +709.5 +609.2 —1316 
139.3 —299.3 *209.8 ¥F44.7 +130.0 248.0 +1698 +246.3 ¥ 404.1 
ci n +26 +25 +33 
D +0.8 +18 —20 
+0.4 +9 ¥6 
u +16 
n —16 —6 —22 
p —25 +31 
+13 
u -19 +3 +13 
v +19 —2 +8 


Norte: The two combinations of relative signs given for each set of coefficients refer to the two sets of sites (e.g., 
Ki,3 and K2,4) related by the twofold screw axis. Experimental errors in the values of the a, }, c coefficients are 
about +0.5 kc/sec for the K sites and about +1 kc/sec for the L and M sites. Experimental errors in the values 
of the n, p, r, u, v coefficients are about +1.4 kc/sec. 
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tensor components the system which are shown Table II. 
The relative signs the off-diagonal components were not uniquely estab- 
lished the coefficients alone because the fact that the rotation 
involved rotation about twofold symmetry axis. the case the sites 
and the relative signs the off-diagonal components could 
determined simply from the signs the coefficients obtained from the Fourier 


TABLE 


co-ordinate system, for the sites colemanite room temperatures 


—38.0 +15.3 +22.8 +210 +139 +299 
+71.8 +25.8 +248 +44.7 130 
1284 +1350 —66.5 +404.1 1698 —246 


Note: The two combinations of relative signs given for each set of tensor components refer to the two sets of 
sites (e.g., and related the twofold screw axis. Experimental errors are about for the 
sites and about +2 kc/sec for the L and M sites. 


analyses the separations the centroids the satellite frequencies from 
the unperturbed frequency. For the and sites, this simple procedure was 
not possible because the second-order effects were small that the shifts 
the centroids the satellites could not detected with the accuracy 
which worked. This problem was solved making few extra measure- 
ments crystal orientations such that none the axes were 
perpendicular parallel the magnetic field. The experimentally determined 
frequencies for these orientations were then compared with the values calcu- 
lated the basis the two alternative combinations signs for the off- 
diagonal components and the correct combination became apparent 
immediately. 

The tensors, were then diagonalized yield the values for the quadrupole 
coupling constants and the asymmetry parameters listed Table III. The 


TABLE 


Quadrupole coupling constants and asymmetry parameters the sites colemanite 
room temperature and —40° 


Quadrupole coupling constant, 


Mc/sec Asymmetry parameter 

Sites Room temperature —40°C Room temperature 

0.566+0.008 

0.862+0.01 
0.825+0.003 

0.359+0.002 

2.553+0.008 
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direction cosines the principal axes the tensors with respect the 
and axes are given Table IV. 


TABLE 


Direction cosines, with respect the and axes, the 
principal axes (x, the electric field gradient tensors the 
sites colemanite room temperature 


Principal axes 


Reference 

+0.247 +0.506 
0.545 +0.574 +0.611 

+0.173 
0.834 


+0.040 +0.550 


Note: Only the relative signs within any given set of the above direction 
cosines are significant. The two combinations of relative signs given for each 
set of cosines refer to the two sets of sites (e.g., Ki,z and Ke,4) related by the two- 
fold screw axis. The errors are about +0.001 in all the values. 


Temperature Dependence the Resonance Lines 

The primary purpose this phase the investigation was find out 
whether n.m.r. techniques are sufficiently sensitive detect phase transitions 
which the atomic rearrangement very small and, so, study the effect 
the transition the electric field gradients. known that some sym- 
metry elements must lost when the crystal transforms. The progress 
transition which symmetry element lost can best examined cole- 
manite orienting the crystal with its axis perpendicular the magnetic 
field and observing the +3/2 +1/2 transitions. Such orientation 
has two advantages. Firstly, can made very accurately, about 
minutes arc, making use our knowledge the symmetry elements 
and the spectrum. Secondly, the nuclei within each group symmetry- 
related sites appear have identical surroundings, and therefore give rise 
single line for each the +3/2 +1/2 transitions long the 
point group symmetry 2/m but soon one more symmetry elements 
are lost each single line will split into two four components (if they are all 
resolved). should possible establish this splitting much more definitely 
than establish slight frequency shift all the lines. 


the temperature varied from about down about 
the boron resonance lines change very small amount frequency but 
remain narrow indicating that the crystal rigorously maintains the point 
group 2/m. our earliest studies, observed that each line starts 
broaden rapidly and about can resolved into two distinct 
components. the temperature further decreased, the frequency separation 
between pair components, which originated from common line, increases 
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until about —14° reaches value about kc/sec (for the site 
149°) after which remains nearly constant until additional 
effect, which will discussed later, occurs. The above sequence events 
appeared perfectly reversible and many thermal cycles such described 
above had measurable effect the room temperature spectrum. The 
separations the satellite lines and originating the and sites 
for the orientation 149° are plotted function temperature 
Fig. Separations between the satellites have been plotted rather than the 
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Fic. Temperature dependences, the region the phase transition, the separations 
between the satellite lines originating the and sites. The crystal orientation was con- 


frequencies the satellites themselves because this emphasizes the change 
due the transition and eliminates certain errors such those due slight 
changes The fact that each satellite line splits into two components 
indicates that the center symmetry has been lost but that the crystal retains 
either the twofold screw axis the glide plane. fact, will explained 
later, colemanite this phase has the symmetry point group 

During our detailed study (described the next section) 
colemanite its phase with point group found that after considerable 
time lapse each the two components (of the original single line) itself 
split into two subcomponents. This, will discussed later, has been taken 
indicate the existence triclinic phase with the point group The 
transformation this triclinic phase was observed occur various tem- 
peratures below about —30°C, the particular temperature which the 
transformation occurs given case depending upon the immediate past 
thermal history the crystal. Once this transformation has taken place, the 
crystal remains its triclinic phase the temperature increased until, 
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-100 -80 -60 -20 


Fic. dependences the separations between the satellite lines originating 
the site. The open circles with the dotted and dashed lines represent the separations 
the crystal transformed first the non-centrosymmetrical monoclinic phase and finally 
the triclinic phase. The dots and solid lines give the separations the crystal transformed 
directly the triclinic phase. All separations were determined constant crystal orientation 
of Oy = 75 . 


the region the four subcomponents, which were originally derived 
from single line, again coalesce. However, the resulting single line some- 
what broader than was before the thermal cycle and this increased breadth 
persists the crystal returned room temperature. With the crystal 
maintained room temperature, the widths the resonance lines gradually 
decrease their original values. the crystal then cooled, the trans- 
formations take place the order described above. However, should 
cooled while the resonance lines are still broadened, the crystal transforms 
about directly from the centrosymmetric the triclinic phase 
which then persists —120° the lowest temperature available us. For 
the sake comparison, Fig. gives reproductions the derivatives 
absorption line originating site, taken first the crystal transformed 
the monoclinic phase and then, after the temperature had been recycled, 
the crystal transformed the triclinic phase. plot the frequency 
separation between the satellite lines originating the sites shown 
function temperature Fig. for transformations both the mono- 
clinic and triclinic phases. Note that when the crystal transforms first 
the monoclinic phase, the later transformation the triclinic phase occurs 


selected temperatures with constant crystal orientation 75°. (a), (c), and 
(d) were recorded during the transformation the non-centrosy mmetrical monoclinic phase 
whereas (a), (f), and (g) were recorded during the transformation the triclinic phase. 
The frequency scale differs slightly from trace trace but about kc/sec per chart division. 
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various temperatures. The splitting the monoclinic phase was studied with 
crystal orientation 149° whereas the orientation 75° was used 
for the investigation the transformation the triclinic phase. 


Figure 9A, reproduction typical spectrum the monoclinic 
phase colemanite —40° shows very clearly, when compared with 


Fic. Selected recordings the derivatives the absorption lines colemanite 

(A) typical spectrum obtained —40° was taken with crystal orientation such 
that the axis made angle 125° with and the and axes ide angles and 
94°, respectively, with the normal the plane containing and the axis. The frequency 
scale approximately kc/sec per chart division. 

(B) example the reduced spectrum —40° obtained the crystal orientation 
150°. The frequency scale approximately 138 kc/sec per chart division. 
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Fig. 2A, that the center symmetry has been lost. Note that all the lines 
are not resolved this particular orientation. stated the previous 
section reduced spectrum again obtained whenever the direction 
parallel perpendicular the magnetic field. Such reduced spectrum 
shown Fig. and this compared with Fig. the equivalent 
spectrum taken room temperature. The frequencies the lines for the 
three rotations have been plotted Figs. 10, 11, and 12. Due the need 
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Fic. dependence angular position the resonance lines colemanite 
—40° for the rotation. 


for larger gap admit the cold cell, less intense magnetic field was used 
this investigation. Thus was 8.985 Mc/sec for the work 
compared with 11.981 Mc/sec for the work room temperature. This change 
which vary inversely were not used except the case the sites 
where only their signs were considered. 

comparison Figs. 10, 11, and with their room temperature counter- 
parts reveals that each room temperature line, with the exception 
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8-800 


Fic. The dependence angular position the resonance lines colemanite 
—40° for the rotation. 


split into two components the crystal transforms. The separations between 
the components derived from and were small that general 
the components could not resolved and appeared two lines shown 
Figs. 10, 11, and except few crystal orientations where the splitting 
could indeed observed. The spectral lines —40° can classified 
exactly the same way the room temperature case. follows from such 
classification that the n.m.r. spectrum colemanite —40° actually 
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contains lines; only these lines are normally observed because eight 
central lines from the and sets are unresolved and appear single 
line and the four central lines from the sites usually appear only 
two lines. conclude that the unit cell colemanite the monoclinic 
phase contains boron sites which the electric field gradients differ 
some respect. Incidentally, the fact that there are boron sites detectable 
the n.m.r. method this phase colemanite compared with the 
sites detectable room temperature confirms our earlier conclusion that 
colemanite centrosymmetric room temperature. 

The quantitative analysis the angular dependence the n.m.r. lines 
has been performed exactly the same way that already described for 
the room temperature case. Since this phase colemanite has point group 
the pairs sites and and are related the twofold axis and 
therefore possess numerically identical VE’s which differ only the orienta- 
tions their principal axes. For this reason, the magnitudes the Fourier 
coefficients for symmetry-related pairs have been averaged and are shown 
Table with their correct relative signs. consequence the fact 


TABLE 
Averaged experimental values the Fourier coefficients (in kc/sec) 
the satellite lines colemanite —40° 


+45.4 —300.2 
+25.0 —310.2 
—43.9 +56.9 54.4 
+106.5 
+141.3 
y —695.4 +614.3 +263 


Note: The two combinations of relative signs given for each set of coefficients 
refer to the two sites (e.g., K: and Kz) related by the twofold screw axis. 
Experimental errors are about +1.5 kc/sec. 


that the rotation was taken about the symmetry axis, again encountered 
the problem assigning the appropriate signs the off-diagonal components 
the case the sites, simple comparison Figs. 10, 11, and 
with their room temperature counterparts (Figs. and shows that 
corresponding off-diagonal components must have the same sign. the case 
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the and sites, such comparison did not give completely unambiguous 
answer one additional measurement was taken order establish the 
signs. The tensor components listed for all the sites Table VI, were 
diagonalized yield the values for the quadrupole coupling constants and 


TABLE 


Average values (in kc/sec) the components the tensor the 
co-ordinate system, for the sites colemanite —40° 


+4.94 +41.4 +204 +116 +300 
+19.1 +15.6 +217 +161 +310 
+86.6 +12.1 —98.7 +199 +54.4 
+54.9 +35.0 —89.9 +299 +32.8 
—1310 +1389 —78.6 1650 263 
1256 +1323 —67.1 +360.0 1742 —242 


Note: The two combinations of relative signs given for each set of tensor components refer to the two sites 
(e.g., Ki and Kg) related by the twofold screw axis. Experimental errors are about +3 kc/sec. 


the asymmetry parameters. These have been listed Table III make con- 
venient comparison with the room temperature results. Table VII gives 
the direction cosines the principal axes the tensors with respect 
the and directions. 


TABLE VII 


Direction cosines, with respect the and axes, the 
principal axes (x. the electric field gradient tensors the 
sites colemanite —40° 


Principal axes 


Reference 
+0.830 +0.289 +0.476 
—0.831 —0.206 —0.517 
—0.566 —0.563 
+0.622 
+0.469 +0.590 +0.657 
—0.544 +0.116 
+0.018 +0.993 +0.113 
+0.559 +0.823 


NOTE: The two combinations of sition signs given for each set of cosines 
refer to the two sites (e.g., Ki and K2) related by the twofold screw axis. The 
errors are about +0.002 in the values for the K and M sites, and about +0.003 
in the values for the L sites. 
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Proton Resonance Line Width 

have carried out preliminary study the temperature dependence 
the proton line width colemanite test the suggestion made Christ, 
Clark, and Evans (1958) that the ferroelectric behavior colemanite can 
traced ordering certain hydrogen bonds which they believe 
state dynamic disorder the room temperature phase. their sug- 
gestion correct, one would expect broadening the proton resonance 
signal the transition temperature. 

The second moment the proton signal plotted function crystal 
temperature Fig. 13. Within our experimental accuracy, the second 
moment did not depend whether the crystal transformed the monoclinic 
the triclinic phase, indicated the signals, the results have all 
been plotted Fig. 13. This plot shows clearly that the proton signal, which 
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Fic. 13. Temperature dependence the second moment the proton resonance line 
colemanite. The experimental error associated with each point about +35 kc?/sec?. 


rather complex, does begin broaden rather abruptly the region 
the temperature further decreased, the proton signal continues 
broaden very slowly until the region about —100° again broadens 
abruptly. There also some suggestion small abrupt broadening the 
signal the region —70° but the experimental accuracy was not good 
enough decide this point. 


DISCUSSION 


The symmetry properties colemanite its different phases have been 
studied other workers using various techniques. The X-ray diffraction 


700 
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study Christ, Clark, and Evans (1958) showed that centrosymmetric 
structure prevails room temperature and they were unable find the 
slightest indication any departure from the space group even 
greatly overexposed films. From study the patterns formed electrical 
breakdown paths, Davisson (1956) concluded that colemanite shows the 
symmetry point group 2/m 150°C and 30°C, and perhaps also 
195° found, however, that colemanite both pyroelectric and piezo- 
electric. His plot the pyroelectric constant against temperature gave 
high and narrow peak —2.5° with appreciable tail the low-tempera- 
ture side and small and rapidly decreasing value for the constant the 
high-temperature side. conclusion was that the essential framework 
the crystal remains centrosymmetrical all temperatures and that the 
asymmetry may due slight displacements certain constituents, 
presumably the light atoms, from centrosymmetrical positions. Chynoweth 
(1957) has also measured the pyroelectric properties colemanite and found 
second-order transition —6°C ferroelectric phase. concluded, 
however, disagreement with Davisson, that temperatures higher than 
about ideal single crystal colemanite would show pyroelectric 
effect and centrosymmetrical above this temperature. 

our own work have found that the unit cell colemanite contains 
non-identical boron sites temperatures below but only distinguish- 
able sites above This indicates quite unambiguous manner that 
the room temperature phase centrosymmetrical agreement with the 
conclusions Christ a/. and Chynoweth. The combination the center 
symmetry with the twofold axis and/or the plane gives the 
point group 2/m which consistent with the space group deter- 
mined X-ray analysis. 

the transition, would expect the space group symmetry reduced 
the basis the structure determination alone. the ferroelectric phase the 
polar axis must parallel the axis the space group 
pendicular the axis the space group but may have any orientation 
the space group Pl. Davisson, Goldsmith (1956), and Chynoweth all 
found the polar axis parallel the axis which rules out the possibility 
Pa. does not entirely rule out the possibility because the departure 
the polar axis from parallelism with the axis may too slight, the 
order few minutes arc, detected the techniques used the 
above workers. Our work demonstrates quite clearly that colemanite may 
transform into phase which only the boron sites are unique and 
this rules out, for this particular phase, the possibility Pl. Thus, 
temperatures below about colemanite has phase which possesses the 
space group This agreement with some recent X-ray work Black 
and Burley, the U.S. National Bureau Standards (reported Christ, 
Clark, and Evans (1958)). 
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reported earlier section, have observed that, with the colemanite 
crystal oriented with its axis perpendicular the magnetic field, single 
resonance line, arising from the same nuclear transition the four symmetry- 
related sites the room temperature phase, may split into two components 
or, under certain conditions, into four components the temperature 
decreased below about The case where splits into two components 
quite unambiguous and has been discussed the previous paragraph. There 
are three possible mechanisms which would lead the single line splitting 
into four components. The first possibility that the crystal has transformed 
into the monoclinic phase with the space group and the same time 
the crystal has become misaligned that the symmetry axis axis) 
longer perpendicular the magnetic field. This would simple instrumental 
difficulty and interest. After numerous tests, all which indicated that 
the crystal had not become misaligned, feel 90% certain that this possibility 
can discarded. The second possibility that the crystal has transformed 
into the phase with space group and the same time twinned. This 
very often happens when ferroelectric crystal transforms. such were the 
case, would expect the four components further split into least 
eight components soon the axis tilted from its position perpendicular 
the magnetic field. The four components were studied various arbitrary 
crystal orientations but additional splitting ever appeared. The remaining 
possibility that the crystal has transformed into heretofore unreported 
phase with the space group Pl. The orientation this new phase not 
uniquely determined the orientation the preceding room temperature 
phase sometimes observed the maximum splitting about 79° 
with minimum about 150° whereas after another temperature recycle 
the minimum splitting would sometimes occur about 79° with the 
maximum about 150°. These angles could not determined very 
precisely the total splitting very small and changes very slowly with 
angle. 

thus appears that the monoclinic phase actually metastable with 
respect triclinic phase. The activation energies associated with these 
transitions would very small indeed but apparently the monoclinic phase 
slightly favored. think very likely that the triclinic phase also 
ferroelectric and that the existence these two alternate phases may offer 
partial explanation for the variety behaviors reported Wieder (1959) 
for different samples colemanite. 

Let now discuss the phase transition. the temperature approaches 
that the transition, selected signal broadens rather slowly which implies 
very small unresolved splitting, then splits into two (or four) components. 
The separation between the components increases fairly rapidly first and 
then approaches value which remains nearly constant the temperature 
further reduced. Such behavior illustrated Figs. and should 
noted particular that the changes the n.m.r. spectrum occurring 
the transition are not due gradual growth the signals arising the 
new phase the expense the signal belonging the centrosymmetric 
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phase. The smooth manner which the lines split the crystal transforms 
suggests that the entire crystal transforms once and that just the 
transition temperature (which cannot, its very nature, exactly deter- 
mined) the atomic arrangement either the low-temperature phases 
virtually identical with that the room temperature phase. the tempera- 
ture lowered below that the transition, the atomic configuration rearranges 
itself until stable arrangement reached which then changes very little 
the temperature further reduced. This behavior very much what 
would expect second-order displacive transition. 

The question remains whether this splitting due actual 
reflection ordering the hydrogen bonds. Since the broadening the 
proton line begins the same temperature the splitting the resonance 
lines, there seems little doubt that the two effects are connected. Because the 
temperature dependence the broadening the proton lines the same 
whether the crystal, indicated the spectra, transforms the tri- 
clinic monoclinic phase and different from the temperature dependence 
the splitting the boron lines, believe that two different mechanisms 
are involved: probably the initial ordering the hydrogen bonds removes 
slightly more stable configuration becomes possible. other words, the order- 
ing the hydrogen bonds may said trigger the reorientation the 
configuration. The broadening the proton line probably 
due the ordering least two groups hydrogen bonds and may 
argued that the splitting the lines reflects ordering one type 
hydrogen bond but that the other type hydrogen bond too far away 
not intimately enough connected cause any effect the boron sites, thus 
explaining the difference between the temperature dependences the proton 
and spectra. According the system hydrogen bonds proposed 
Christ ordering possible for the hydroxyl group and likely for the 
water molecules, Os. Any ordering the group would expected affect 
the boron atoms most strongly since forms part the 
chain. connected the chain only through hydrogen 
bonding one the oxygen atoms, O;. Thus were miss the effects 
ordering process among these hydrogen bonds would expect 
the one involving the water molecules; yet Christ al. feel most strongly 
that this the very ordering process which leads the ferroelectric behavior 
and therefore must occur about where observe the major effect 
the lines. possible that thorough study the complex proton 
resonance will help clarify this point. 

principle, should possible use the quantitative data the 
quadrupole coupling tensors the room temperature and ferroelectric phases, 
derive the structural changes which take place the transition. However, 
the computational procedures are likely prohibitively difficult since 
this case, one cannot use simple model point charges. likely that 
information about the structural changes can best obtained neutron 
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diffraction study colemanite although additional difficulty, that the 
high-absorption cross section the nuclei for thermal neutrons, will 
encountered. 

The transition temperature colemanite has been variously reported 
(Goldsmith 1956), (Chynoweth 1957), and —7° (Wieder 
1959). have found that our crystal the transition occurs temperature 
several degrees higher than any these values but there probably little 
significance this variation. Wieder reported finding wide variation 
transition temperatures, between —7° and for number samples, 
which attributed internal space changes caused impurities. Samples 
from the same collection colemanite from which our crystal was taken 
give biased hysteresis loops which suggest that the transition our crystal 
may have been elevated space charge effects. One interesting point regarding 
the transition temperature the fact that observed the transition the 
triclinic phase take place consistently higher temperature, than 
the transition the monoclinic phase This suggests that two separate 
transitions are involved which fortuitously occur very close together and 
the room temperature phase happens undercooled below the temperature 
the first transition transforms directly the monoclinic phase which 
metastable. can course argued that cannot define exact transition 
temperature because the very nature the transition and that simply 
observe the change the triclinic phase higher temperature because the 
final splitting larger than the monoclinic case. Following this argument, 
there only one transition with two alternatives, but only slightly different, 
boron-oxygen configurations available the structure. The first broadening 
the proton line not abrupt enough for able set sufficiently 
small error the transition temperature help distinguish between the two 
alternatives and have been unable settle the point. 

certain extent, one may correlate the boron sites, determined 
the n.m.r. results, with the structural co-ordination the boron atoms. 
the discussion follow will assumed that only nearest neighbors make 
significant contributions boron site. Christ al. have shown that 
two the three boron atoms the asymmetric unit (at room temperature) 
are located the interstitial positions near the centers two slightly dis- 
torted oxygen tetrahedra and that the third occupies the interstitial position 
the center slightly distorted equilateral triangle oxygen atoms. 
the center perfect tetrahedron the symmetry cubic that must 
zero. slightly distorted tetrahedron, however, one would expect the 
magnitude and therefore the quadrupole coupling constant small 
the center although need longer zero. Without detailed knowledge 
the distortion and resulting charge density distribution one cannot predict 
the orientation the principal axes nor the value the asymmetry 
parameter. The symmetry the center perfect equilateral triangle 
requires that zero and that the principal axis parallel the normal 
the plane the triangle but puts restriction the magnitude VE. 
Thus the case the BO; triangle, would expect the quadrupole coupling 
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constant considerably larger than for the BO, tetrahedra and, unless 
the triangle greatly distorted, that the value should small and 
that the principal axis should nearly perpendicular the plane con- 
taining the three oxygen atoms. shown Table III, have indeed 
found that the quadrupole coupling constants are small the and boron 
sites, whereas the quadrupole coupling constant considerably larger and 
very small the site. Also, follows from the atomic co-ordinates 
(Christ 1958) that the principal axis the electric quadrupole coupling 
tensor the site nearly parallel the normal the plane containing 
the oxygens Os, and Fig. the triangle. This clearly shown 
the direction cosines with respect the and directions, 
respectively, shown below: 
direction cosines principal axis 


direction cosines normal plane 
+0-109 


There seems little doubt that the boron sites which have labelled and 
are the oxygen tetrahedra and that the site the oxygen triangle. 
have been unable find unambiguous correlation between our and 
sites and the individual boron sites the tetrahedra. 

Incidentally, interest note that the n.m.r. study the boron 
sites kernite (Waterman and Volkoff 1955) revealed that the electric field 
gradients existing sites and (using Waterman and nomen- 
clature) are very similar those our sites and while the electric 
fields existing boron sites and are very similar those our site 


TABLE VIII 


comparison the values obtained Waterman and Volkoff for the quadrupole coupling 
constants and asymmetry parameters for the sites kernite with the values obtained 
for colemanite room temperature 


Quadrupole cou- 645 588 2563 2567 Quadrupole cou- 436 
pling constant, pling constant, 

Asymmetry para- 0.54 0.60 0.163 para- 0.487 0.825 0.059 


meter meter 


comparison given Table view this striking similarity, suggest 
that kernite contains two BO, tetrahedra and two triangles, contrary 
the X-ray analysis results Portoles (1948). This suggestion has been sup- 
ported Dr. Christ (private communication), who has pointed out that the 
chains should the same kernite borax. Morimoto 
(1956) has shown that borax, which has the same chemical formula kernite 
except for the water hydration, contains two BO, tetrahedra and two 
triangles its asymmetric unit. 


A 
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RELATIVE ENERGIES TILT-TYPE 
SUBBOUNDARIES ALUMINUM! 


Aust 


ABSTRACT 


The relative energies subboundaries were obtained, terms reference 
large-angle grain boundary, aluminum bicrystals having two different 
effective impurity concentrations. The data for subboundary energy versus 
subgrain orientation difference (@) showed cusp 0.75°. The existence 
such cusp was previously predicted Martius and Chalmers from con- 
siderations impurity effects boundaries. 


INTRODUCTION 

During study grain boundary migration motivated subboundaries 
bicrystals aluminum (Aust 1959), became apparent that the relative 
energies the subboundaries could obtained terms the grain boundary. 
The method consisted measuring the dihedral grain boundary angle the 
point intersection each subboundary with the grain boundary, after 
annealing near the melting point the metal. The results, described 
this report, indicate the presence energy cusp orientation difference 
about 0.75°. The existence such cusp was previously predicted 
Martius and Chalmers (1952) the basis impurity effects boundaries. 
This observation believed represent the first experimentally measured 
energy cusp for interfaces, except for those existing twin boundaries. Several 
the consequences this phenomenon, relation properties boundaries 
metals, are discussed. 


EXPERIMENTAL PROCEDURE AND OBSERVATIONS 

Single crystals zone-refined aluminum, approximately sq. and 
length, were obtained zone-melting aluminum 99.996 wt.% starting 
purity graphite boat under helium atmosphere. molten zone 
length, which was maintained induction heating, was passed along the 
length the bar usually rate 0.9 per minute, except the last pass 
which the rate was 0.2 per minute. Two single-crystal specimens having 
different impurity contents were produced passing the molten zone along 
the entire length the bar for passes one case and passes another. 

Direct measurements the 4-pass and 12-pass aluminum single crystals 
using the method Bean, DeBlois, and Nesbitt (1959) gave the following data 
for the resistance ratio between the electrical resistance room temperature 
and the resistance liquid helium temperature: 

4-pass crystal, center bar 3100, 
12-pass crystal, center bar 9000. 

received November 17, 1959. 

Contribution from the General Electric Research Laboratory, Schenectady, New York. 
Can. Phys. Vol. (1960) 
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The resistivity ratio for the center portion the 12-pass aluminum may 
the highest yet reported for aluminum, but still represents effective impurity 
content few parts per million. addition, the thermally activated motion 
grain boundaries, which known very sensitive solute impurity 
content (Rutter and Aust 1958), was also studied. was found that the 
velocity grain boundaries was orders magnitude faster, and with 
much smaller temperature dependence, the case the 12-pass aluminum 
compared with the 4-pass aluminum (Aust 1959). 

evident from both the electrical resistivity and grain boundary migration 
experiments that the effective impurity content the 12-pass aluminum 
crystal less than that the 4-pass crystal. 

Single crystals about length were removed from the center portion 
each zone-refined bar cutting them with jeweller’s saw, then etching 
them aqua regia remove any strained material the cut surfaces. The 
single crystals were electrolytically polished 2:1 solution methyl alcohol 
and concentrated nitric acid, using intermittent polishing periods for total 
time minutes current density amp per sq. dm. The 
specimens were then chemically polished Alcoa R-5 solution. X-Ray and 
micrographic studies indicated the presence striations lineage sub- 
structure the single crystals, having misorientations high 2°. The 
orientations the 4-pass and 12-pass aluminum single crystals, containing 


substructure, are given Fig. 


oT “00! 


Fic. Stereographic plot specimen axis growth direction and normal top 
surface starting single crystals; 4-pass (light) and 12-pass (shaded). 


Bicrystal specimens, containing one crystal with lineage and the other 
crystal without lineage, were next obtained from the single crystals using the 
following technique, previously described Rutter and Aust (1958, 1959). 
One end the striated single crystal was deformed locally compression and 
the specimen then annealed argon Numerous new recrystal- 
lized grains are thereby introduced one end the specimen, from which 
sometimes found that one recrystallized grain grows preferentially into 
the lineage substructure the undeformed part the crystal. this way, 
was possible obtain single grain produced recrystallization growing into 
striated melt-grown crystal. Figure shows example such growth the 


12-pass aluminum. 


PLATE 


Fic. Photomicrograph showing grain boundary subboundary interaction. 1000X. 
Fic. Macrophotograph showing anchored grain boundary 4-pass bicrystal. 


‘ 
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The orientation differences the subboundaries, such depicted Fig. 
were obtained accuracy +.05° positioning X-ray beam the 
subboundary and thereby recording the Laue reflections from two adjacent 
subgrains the same film. Analysis the X-ray Laue back-reflection films 
revealed that the subboundaries, which run nearly horizontally (Fig. and 
which intersect the grain boundary, were tilt-type boundaries having mis- 
orientations 0.1 2°. Several the subboundaries which are approximately 
vertical (Fig. were also analyzed and found twist-type boundaries, 
with misorientations the order 0.5°. 

The bicrystals were chemically polished prior each argon anneal, and the 
interfaces were clearly revealed thermal etching seen Fig. Each 
triple point, where tilt-type subboundary intersected the grain boundary, 
was photographed X500 X1000 after several annealing treatments 
640° argon. The dihedral angle, made the grain boundary each 
point intersection with subboundary, was measured within +0.5°, 
directly from the negatives from prints enlarged times the original 
magnification. typical example grain boundary interaction 
shown Fig. The relative energies the tilt-type subboundaries were 
obtained terms the reference grain boundary, accuracy +.01, 
using the relationship 


where the energy the subboundary, the energy the reference 
grain boundary, and the dihedral angle opposite the subboundary. 
measurement then yields the value ys, when assumed equal 
unity. The angle varied from 165° 179.5°, resulting values from 
0.26 0.01. 

The relative energy the subboundaries (ys) function the subgrain 
orientation difference was obtained for two bicrystal specimens, one from 
the 4-pass aluminum and the other from the 12-pass aluminum. The orientation 
differences the two reference large-angle grain boundaries are given 
Fig. where the single axis relating the starting striated crystal with the 
adjacent recrystallized grain the smallest amount rotation plotted. 


(DEGREES) 


Orientation relationships for the two bicrystals. 4-pass Al, axis 12-pass Al, axis 


RESULTS AND DISCUSSION 


The results obtained for the relative energy the subboundaries (ys) 
terms large-angle grain boundary, function orientation difference 
the subgrains, are shown Figs. and for the 4-pass and 12-pass 
aluminum bicrystals respectively. The two curves are similar that 


550 CANADIAN JOURNAL PHYSICS. VOL. 38, 1960 
| 


w 
<= 
=> 
> 
= 
| 
w 
2 
0 


DIFFERENCE ORIENTATION (DEG) 


Relative subboundary energy versus orientation difference the subgrains 


4-pass Al. 
Relative subboundary energy versus orientation difference the subgrains 12-pass 


Al. 


rises rapidly with increasing reaching maximum value 
0.7-0.8°, and then increases again with further increase The main 
difference the two curves Figs. and the relative drop from 0.16 
0.06 for the 4-pass aluminum, and from 0.15 0.09 for the 12-pass aluminum, 
i.e., aslightly more pronounced energy cusp present the aluminum bicrystai 


lower purity. 
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The results depicted Figs. and were obtained under essentially dynamic 
conditions where the subboundaries provided the driving force for the motion 
the grain boundary. order determine whether this data represented 
equilibrium values (where equation (1) only applicable), the following 
experiment was conducted. Grooves were acid-cut opposite sides the 
4-pass aluminum bicrystal, just ahead the moving grain boundary. The 
specimen was then annealed 640° argon until the grain boundary became 
anchored the grooves, reaching equilibrium position which the lineage 
driving energy was balanced equal and opposite driving energy associated 
with the energy and curvature the grain boundary. Figure shows the 
anchored grain boundary equilibrium condition after annealing 640° for 
hours. Two the grain boundary dihedral angles, points intersection 
with the subboundaries, were then remeasured and the relative energy the 
subboundaries obtained from equation (1) before. The results for the 4-pass 
specimen were follows: 


(deg) 
subboundary Dynamic Equilibrium 
0.06 0.07 
0.16 0.19 


These results indicate that the energy cusp, found under conditions 
moving grain boundary, also present under equilibrium conditions. 
previous study grain boundary energies tin tricrystals (Aust and Chalmers 
1950), the size the dihedral angle opposite the grain boundary, whose 
orientation difference was varied from specimen specimen, was studied 
function annealing time addition, the rate movement 
the triple point the intersection three grain boundaries was also deter- 
mined. was found that the angles between the grain boundaries reached 
equilibrium values well before the triple point reached its equilibrium position. 
This observation also indicates that grain boundary angles measured from 
grain boundaries which can move further annealing may represent equili- 
brium values the triple points. 

The boundary energy cusp shown Figs. and was predicted Martius 
and Chalmers (1952), based considerations the free energy system 
which contains impurity atoms and boundary. Their preliminary calculations 
for the possible dependence relative boundary free energy upon difference 
showed minimum the energy about 0.5°; however, 
the exact calculated value depends the assumption the magnitude 
the volume strain field the individual dislocations boundary. These 
calculations are based dislocation model using elasticity theory for the 
strain around the dislocation and the model the core 
distortion (Martius 1959). The calculated value consistent with the 
observed value about 0.75° and suggests that the same phenomenon 
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responsible both cases. The energy cusp would expected absent 
the case pure metal with impurity the boundary. The present 
experimental data indicate decrease the magnitude the energy cusp 
with increasing purity the aluminum, suggesting that the observed energy 
cusp may related impurity effects considered Martius and Chalmers 
(1952). Since impurities apparently exert greater influence the sub- 
boundary energies the cusp position 0.75°, the energy the latter 
value should then the one most sensitive varying impurity content. 
This may account for the fact that the main difference the curves Figs. 
and for the 4-pass and 12-pass aluminum respectively, occurred the 
cusp position. was pointed out Chalmers (1952), this cusp may 
important factor the persistence lineage substructure metal crystals 
which are grown from the melt. The formation and/or retention such 
boundaries crystals would energetically more favorable when impurities 
are present. 

recent study absolute interfacial energies copper having purity 
99.98% showed evidence cusp the boundary energy versus orientation 
difference for lineage boundaries with misorientations from 0.65 
(Gjostein and Rhines 1959). However, possible that cusp, present 
copper, may have been missed since the values are given +0.2°, that 
cusp exists value below 0.65° this relatively more impure material. 
One interesting consequence the energy cusp found aluminum that 
subboundary with value about 1.6° could form two subboundaries, each 
with value 0.8°, but with total interfacial energy than the single, 
higher-angle subboundary. energetic basis, therefore, subboundaries 
with between 1.4 and 1.8° may then absent; this was found the case 
the present aluminum bicrystals (e.g. see Figs. and 7). 

The present energy cusp was not predicted the theoretical treatment 
Shockley and Read (1950) for dislocation boundaries, which the effect 
impurities was essentially neglected. The results Figs. and are given 
single plot relative boundary energy versus Fig. addition, two 
curves are shown which were calculated from the Shockley and Read equation 
that the exact position the two calculated curves depends the magnitude 
which the boundary energy assumed unity. The latter value 
used here was 20° although this value has very little effect the curves 
small angles 

The Shockley and Read curve with similar that frequently 
observed relative grain boundary energy experiments impure materials. 
However, the curve with —0.5 describes the present data the region 
about 0.3° which the boundary energy least influenced 
impurity due the larger spacing the dislocations the boundary (Chalmers 


1959). The curve with —0.5 might expected, therefore, represent 
better the simple dislocation model the boundary, without impurities, than 
the curve with +0.5. fact, negative value for —0.3 may 


obtained from the Shockley and Read calculations (Brooks 1952). addition, 
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Fic. Composite plot relative energy versus data for the subboundaries the 4-pass 
and 12-pass Al, together with theoretical curves (dashed) with —0.5 and 
+0.5. 


the magnitude the orientation difference, where the relative boundary 
energy reaches maximum, may calculated from The @nax 
values are then 13° and 35° for —0.5 and +0.5 respectively. The dis- 
location model the boundary probably suitable description the 
boundary structure for values approaching 13°, but not for values high 
35°. 

The above discussion Fig. indicates possible explanation the previous 
agreement the experimental data grain boundary energies with the 
Shockley and Read calculations far beyond the range where the theory 
longer applicable. This explanation based the suggestion that the boundary 
energy metal may rise more rapidly with increasing than that 
metal. the case the usual relative boundary energy measure- 
ments tricrystals, this argument assumes that the energy the reference 
large-angle boundary less influenced impurity than that the low- 
medium-angle boundary whose relative energy measured. This assumption 
appears reasonable since the energy low-angle and medium-angle 
boundaries (i.e., 15°) markedly dependent structure, while the 
energy large-angle boundaries (i.e., 15°) independent structure, 
except for twin boundaries and other possible special orientation relationships. 

evident from the present experiments subboundary energies, and 
from those grain boundary migration (Rutter and Aust 1958, 1959), that 
the role impurities the properties boundaries should not neglected. 


SUMMARY 


The relative energy subboundaries was studied function orientation 
difference from 0.1 aluminum bicrystals having two different effective 
impurity contents. cusp was found the energy versus relationship 
0.75°, which agreement with the predictions Martius and Chalmers 
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based impurity effects boundaries. The significance the cusp with 
regard substructures metal crystals and previous work boundary 


properties discussed. 


ACKNOWLEDGMENTS 


The zone-refined aluminum used this study was produced collaboration 
with Hilliard. The electrical resistivity data were obtained 
DeBlois and Felton. The author would also like thank Chalmers 
for his helpful comments the manuscript this report. 


REFERENCES 


Aust, and 1950. Proc. Roy. Soc. 201, 210. 

Aust, and Rutter, 1959. Trans. Am. Inst. Met. Engrs. AIME 215, 119. 

P., DEBLots, W., and 1959. Appl. Phys. (12), 1976. 

Brooks, 1952. Metal interfaces (American Society for Metals, Cleveland), 20. 

1959. Private communication. 

1959. Private communication. 

and 1952. Discussion imperfections nearly perfect crystals 
(John Wiley Sons), 371. 

Reap, and 1950. Phys. Rev. 78, 275. 

and Aust, 1958. Acta Met. 375. 


NEUTRON SOURCES FROM THE BERYLLIUM 
REDUCTION PLUTONIUM 


ABSTRACT 


Sintered neutron sources have been prepared reacting plutonium dioxide 
with beryllium metal powder. The reaction, forming alloy, 
proceeded very slowly 750° but was nearly complete 850° Sintered, 
non-friable agglomerates were obtained heating the reacted mixtures 1250 
1450° depending the beryllium concentration. the composition range 
particular interest, i.e. Be/Pu atom ratios 13:1 and above, the neutron output 
from the sintered products were from 98% the theoretical neutron yield. 
The presence BeO the sinters was found have negligible effect the 
neutron output the alloy. Studies some factors affecting neutron output 
sinter quality are also reported. 


INTRODUCTION 


alloy sources have been prepared the past two 
methods: reducing plutonium fluoride with excess beryllium (Runnalls 
and Boucher 1956) and heating the two metals together (Coffinberry and 
Tate 1958). 

However, the direct reduction with could used, simpler, 
more economic method preparation would result since there would 
need for fluorination metal preparation. 

According published thermodynamic data (Coughlin 1954; Holley 
1958), BeO more stable than 1000° since for the reaction 


the free energy change 28+1 kcal. 
excess present, the equilibrium shifts further the right because 
alloy formation, and the reaction takes the form 
This paper describes procedures that have been developed for the preparation 
neutron sources this method. Also reported are: 
(a) brief study the effect temperature the rate which reaction 
(2) proceeds vacuum. 
The effect plutonium dispersal. 
(c) The effect compacting the reactants before sintering. 
(d) The effect partial melting the alloy. 


(e) The variation the neutron yield the sinters with composition. 
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EXPERIMENTAL PROCEDURE 


(a) Materials 
Plutonium dioxide 99.9% purity was prepared precipitating plutonium 


oxalate from nitrate solution with oxalic acid crystals, then converting the 
oxalate the dioxide heating air 650° for hours. This material 
was found correspond the formula from the loss weight after 
igniting sample constant weight 1300° (Drummond and Welsh 1957). 
The isotopic composition the plutonium determined mass spectro- 
meter analysis was: 
93.558+0.012 atom 
atom 
0.494+0.008 atom 


Results obtained from multichannel pulse analyzer showed that 2.6% the 
total alphas were from which corresponds 0.011 atom Hence, 
the calculated specific activity was d.p.m. per Pu. 

The plutonium dioxide and the 99.8% pure beryllium metal powder were 


both sized —200 mesh. 


(b) Preparation the Samples 

The powdered reagents were placed stoppered bottle and shaken 
mechanical vibrator. For early samples, the operation was carried out near 
neutron counter indicate the increase neutron emission mixing pro- 
ceeded. The neutron output increased rapidly during the first minute but 
observable increase resulted from mixing for longer periods. mixing time 
minutes was considered adequate, therefore, for all subsequent samples. 

The loose mixture was placed recrystallized-alumina crucible, 9/16 in. 
was then degassed heating slowly under vacuum tungsten- 
element furnace 400° hours and was held that temperature for 
minutes. The furnace was then cooled and cylindrical 30-g tantalum slug 
1/2-in. diameter alumina disk was placed top the charge. Compression 
the charge, which had avoided degassing prevent ejection 
powder, was found necessary during sintering obtain uniform sinter from 
top bottom. 

Sintering was carried out heating for minutes vacuum either 
1250° 1450° depending the beryllium concentration, and holding 
the charge this temperature for minutes. 

The choice the sintering temperature was influenced two considera- 
tions: the production satisfactory sinter and the prevention partial 
liquation the alloy. 

satisfactory sinter was defined one that was dust-free, non-friable, and 
strong enough not break normal handling. According the amount 
excess that was used reaction (2), the end products varied from highly 
ceramic highly metallic. For samples with high concentration BeO 
the end products (Be/Pu atom ratios around 15:1) temperature from 1400 
1450° was found necessary produce acceptable sinters. the other 
hand, with highly metallic products (Be/Pu atom ratios 30:1 and upward), the 
desired sinter qualities could obtained heating 1250° 
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was also observed that samples with Be/Pu atom ratios above 30:1 could 
not heated 1450° without causing partial liquation the alloy accom- 
panied decrease 10% the neutron output. This effect, discussed 
more fully later, was avoided sintering samples that range compositions 
1250° 

Since the evaporation loss beryllium vacuum was appreciable 1450° 
the sintering time had kept short. the sintering period was longer than 
minutes, was difficult obtain predictable Be/Pu atom ratio the 
final product. All samples were, therefore, held for minutes the maximum 
temperature. The evaporation losses this period were 1450° 
and less than 1250° 

establishing the Be/Pu atom ratio the sintered samples, the content 
was determined from the weight oxide put the charge. The content, 
which included alloy and BeO, was determined subtracting the 
evaporation loss from the amount weighed for the charge. 


(c) Neutron Counting Method 

For neutron output measurement, the sintered sources were transferred from 
the glove-box into screw-top bottles. The neutron yields were determined 
area low neutron background direct comparison room temperature 
the neutron output each source with that standard Pu-alpha-Be 
source using the simple paraffin-block counter described Hanna and 
Runnalls (1956). 

The standard source had been calibrated the National Research 
Council, Ottawa (Geiger 1958), against the N.R.C. Ra-alpha-Be standard 
(N-200-1) which turn had been compared with the N.B.S. primary standard. 
The neutron output the standard was the same order magnitude 
that the sintered sources made this work, i.e. n/sec. 


’ 


RESULTS 


(a) Effect Temperature the Rate Reduction Plutonium Dioxide 
Beryllium Vacuum 

was determined from X-ray diffraction analysis that practically 
reaction occurred between plutonium dioxide and excess beryllium after heating 
the reactants 750° for minutes vacuum. the other hand, after the 
same time 850° reduction the was found well over 
complete. 

The X-ray diffraction patterns residues from leaches sources 
sintered 1250 1450° showed only the lines BeO, whereas the diffrac- 
tion pattern unsintered sample containing 0.5 and 99.5 
mol.% BeO showed all the lines but one. This would indicate that 
the sinters reduction had proceeded virtually completion. 


(b) The Effect Plutonium Dispersal 

With sources made the conventional methods, the ideal distribu- 
tion the plutonium may approached the constituents are completely 
melted single liquid and cooled rapidly. neutron sources, 
however, attempts melt the alloy resulted separation Be-rich material 
from the sinter and decrease the neutron yield. 
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The sintering operation itself gave only limited increase the neutron 
output, and was therefore desirable achieve mechanically high degree 
dispersal the sample before sintering it. This could possibly best 
done reducing the maximum particle size. 

Experiments were carried out determine the effect various mesh sizes 
the neutron yield after sintering. The samples, all with Be/Pu 
atom ratio 100:1, were sintered 1250° according the method described 
earlier. The neutron yields were determined after cooling with results shown 
Table The dusting loss during the sizing operation, the time taken 


TABLE 
Effect mesh size the neutron yield 


mesh size theoretical 
—100 
200 
—300 


produce —300-mesh were both appreciably higher than for —200-mesh 
material. Since only slight improvement the neutron yield resulted when 
using —300-mesh material, the results obtained with —200-mesh were 
considered satisfactory and this size was adopted for subsequent experiments. 

Sintered samples from the sizing experiments were cross-sectioned and 
polished for microexamination. Because the porosity the samples, was 
found preferable examine them the as-polished condition, without etching. 

The effect mesh size the porosity the samples sintered 
1250° shown low magnification Fig. 1A. The porosity, quite appre- 
ciable the —100-mesh samples, decreased with decreasing coarseness 
the particles and was relatively very low the —300-mesh sample. 

The effect mesh size the size the larger particles pro- 
duced shown high magnification Fig. 1B. The phase appears 
more less circular islands darker material the bright beryllium 
matrix. The irregular black areas are pits caused partly the porosity the 
sample and partly the removal material polishing. The diameter the 
largest particles observed was about microns the —100-mesh 
sample, microns the —200-mesh sample, and about microns the 
sample. There were very few large and intermediate size particles 
most were very fine. 


(c) The Effect Pressing the Mixed Reagents Before Sintering 

The constituents for nine samples with Be/Pu atom ratios from 13:1 
120:1 were mixed and then pressed 20,000 p.s.i. The resulting green pellets 

Fic. As-polished sections samples 100:1 Be/Pu atom ratio sintered 1250° 
mesh size indicated circles. Low magnification (A) shows the effect mesh 
size porosity. High magnification (B) shows the effect mesh size the size the 
largest particles produced. The light phase matrix, the darker more less 
circular phase and irregular black areas are voids. 


after sintering 1250° for minutes 
Neutron yield, 


@ 


PLATE 


As-polished section sample 32:1 Be/Pu atom ratio which partial liquation 
had started occur while sintering 1450° Melted material (top right) had exuded from 
the sinter. 
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were dust-free and could withstand normal handling without breaking 
chipping. The pressed pellets could heated vacuum minutes the 
sintering temperature without cracking ejection powder. Degassing, 
reduction, and sintering were, therefore, carried out continuous heating 
operation, eliminating the separate degassing step necessary loose powder 
samples. The samples were held for minutes the sintering temperature, 
1250 1450° depending the composition. 

Pressed samples were sintered with and without the tantalum slug weight 
top the pellet. both cases the pellets showed dimensional increases due 
chemical reaction. the presence the slug pressure, the increase height 
and diameter after sintering was less than 2%, while without the pressure, the 
dimensional increases varied from with visible distortion the 
shape. 

The neutron yields the compacted samples after sintering were essentially 
the same those from corresponding unpressed samples. 

Compacted samples have also been sintered tank argon following the 
same procedure with neutron yields generally lower than those 
equivalent sources sintered vacuum. 


(d) Effect Partial Melting the Alloy 

detectable amount melting occurred when samples Be/Pu with 
atom ratios between 7:1 and 30:1 were heated 1450° and held there for 
minutes. However, samples with ratios above 30:1 were heated 1450° 
partial melting would occur and drop 10% the neutron output 
would result. appeared interest determine the reason for this drop 
the neutron yield. 

The polished cross section sample which melting had just started 
take place shown low magnification Fig. This sample had Be/Pu 
atom ratio 32:1 and had been sintered 1450° for minutes. The bright 
Be-rich phase appearing top right melted and separated from the sinter and 
had started flow down into the space between the shrunken sinter and the 
crucible wall. This amount alloy separation was enough cause drop 
the neutron yield. 

order study partial melting over the duration the sintering operation, 
three samples with Be/Pu atom ratio 100:1 were degassed, heated 1350° 
hours, and then cooled. The neutron output during heating and cooling 
was measured neutron counter positioned near the furnace. The neutron 
vields were later related the absolute yields comparison the cooled 
sample with standard source. The time which melting and freezing occurred 
was determined from the change slope the recorded output thermo- 
couple placed near the sample. 

all three samples, two separate approximately equal decreases the 
neutron output were observed. The first started when melting began 1284° 
and continued until was reached where the power was reduced 
allow cooling. The second decrease the neutron yield started with the 
beginning solidification about 1284° and continued until all the melted 
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Variation the neutron emission rate from 100:1 Be/Pu atom ratio sample during 
heating 1350° and cooling room temperature. Partial melting started 1284° 


material had solidified. The neutron output then remained constant room 
temperature. The decrease the neutron yield was approximately during 
melting and the same during solidification. These observations are shown 
graphically Fig. 

One the samples was cross-sectioned and polished and shown low 
magnification Fig. appreciable amount material had melted 
and flowed out the sinter form ring (torus) around the porous center 
portion. The two parts had about the same volume and analysis indicated 
Be/Pu atom ratios 700:1 and 50:1 for the melted and unmelted parts 
respectively. 

Thus, when partial melting occurs, Be-rich phase that becomes liquid. 
The distribution longer uniform and this accounts for the first 
the two observed drops neutron yield. 

The bottom left corner the sample shown Fig. appears five times 
larger Fig. 4(6). this magnification considerable segregation the acicular 
phase can observed the area left indicated arrows. Higher 
magnification portion this area segregation (Fig. 4(c)) shows numerous 
particles the matrix and these have grown considerable size. 
Growth large particles the melted portion, probably resulting 
from slow cooling, could account for the drop the neutron yield that occurred 
during solidification. 

was decided that samples with Be/Pu atom ratios greater than 30:1 
would sintered 1250° This avoided melting, gave good sinters, and 
kept neutron output better than theoretical. 


PLATE 


Fic. As-polished section 100:1 Be/Pu atom ratio sample heated 1350° produce 
partial 
(a) Low magnification shows center portion that remained solid and two bright areas 


Be-rich material which had melted and formed ring around the sample. 

(b) Shows the bottom left corner sample medium magnification where 
the appears the melted material (arrows). 

High magnification the segregation area shows large particles dark the 


light matrix. 


zregation 


fi 
= 
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(e) the Neutron Yield with Atom Ratio 
The neutron yields from sintered sources have been plotted Fig. 
show the effect the Be/Pu atom ratio the observed yield neutrons. 
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Fic. variation the neutron yield with atom ratio for sintered sources made the 
reduction with Be. 


The curve Fig. represents the theoretical neutron yield and was con- 
structed from the expression used Runnalls and Boucher (1956), i.e. 


Mnax = N XSpe/(N X Spe +Spu) 


where fraction the maximum theoretical neutron yield, 
Be/Pu atom ratio, 
Spe atomic stopping power 0.63, 
Spy atomic stopping power 4.6. 


The value the thick target yield, determined Runnalls and 
Boucher (1956) was alphas, when related the Chalk River Ra- 
alpha-Be standard source. mentioned Hanna and (1956), 
appeared that the Chalk River standard, which had been calibrated against 
the Argonne National Laboratory source No. 38, was too high with respect 
the National Bureau Standards Ra-gamma-Be source No. Since 
now generally recognized that the N.B.S. standards are the more accurate, 
the quoted neutron output for the Chalk River standard has been lowered 
recently 5%. Hence, the thick target neutron yield for alphas 
beryllium, has been decreased alphas. 

will noted from Fig. that good agreement obtained between the 
highest observed values obtained this work and the theoretical neutron 
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yields. The neutron yield obtained for mixture produced according the 


reaction 


36.6 alphas which good agreement with the Runnalls and Boucher 
(1956) adjusted value 36.7 alphas for This would suggest 
that the BeO not contributing significantly the neutron output from the 
alloy. 

Practical applications neutron sources very often require minimum 
physical size for given neutron yield. For fixed amount plutonium, 
sources with high Be/Pu atom ratios produce higher neutron yields than those 
with low ratios, but they also tend become too bulky for many applications. 

Konobeevsky (1955) reported that the maximum neutron yield per unit 
weight alloy was obtained Be/Pu composition 13:1, while results 
Los Alamos indicated that higher yields per unit weight alloy were obtainable 
lower content than the composition (Coffinberry and Waldron 
1956). report, however, has been found the literature the variation 
with composition the neutron yield per unit volume alloy. The maximum 
neutron yield per unit volume alloy may expected near the 
composition provided the theoretical densities the alloys are obtained. 
order approach theoretical density the region composition, 
however, temperature 1600° more required. The sinters this work 
were heated only 1450° and was expected, therefore, that the 
density sinters near the composition would appreciably lower 
than the theoretical density. 

became desirable, therefore, determine what composition the maxi- 
mum neutron yield per unit volume sinter was obtained when sintering 
1450° The sinters were coated with wax and their volumes measured 
pycnometric method. The neutron yields Fig. were converted neutrons 
per second per cubic centimeter sinter and results were shown Table 


TABLE 


Variation the neutron yield per unit volume 
products sintered 1450° for minutes 


Neutron yield, 


Density, n/sec 
Be/Pu atom ratio sinter 
9.3 1.624 1.03 
11.0 1.940 1.03 
12.6 1.870 1.05 
1.741 0.98 
14.2 1.660 
15.0 1.680 0.95 
16.7 1.901 1.07 
17.4 2.019 
20.5 2.240 1.28 
2.141 
27.8 2.189 
32.0 2.229 1.08 
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maximum value indicated Be/Pu atom ratio the region 20:1. 
Due marked difference the sinterability the various compositions 
1450° samples 15:1 atom ratio reached only 40% the theoretical density 
while samples 20:1 ratio reached 60%. 

DISCUSSION 

From series more than measurements the yields two standard 
sources, the possible errors the comparative neutron yield measurements 
were estimated +2.5%. This subject further deviation 
include the error absolute standardization and intercalibration measure- 
ments the primary standard. 

The stability neutron sources relies the long half-life 
years). However, the isotopes and Pu*! are always 
present varying amounts with depending the characteristics the 
reactor and the length time the plutonium was irradiated. Since the alpha 
half-lives and are fairly long, years and 6600 years respectively, 
the decrease the neutron output from decay these isotopes slight; for 
the piutonium used this work, was calculated that the decrease the 
alpha emission rate would less than 0.5% years. 

however, decays beta-particle emission with half-life 
the alpha emitter and present appreciable amounts, 
gradual increase the neutron output from sources will result. Since 
the plutonium used this work was rather high growth correction 
factor was required and, calculated Hanna (1958), amounted 1.4°% 
per year. The sources prepared reduction have been found 
satisfactory for many applications such instrument calibration, low 
flux experimental work A.E.C.L. special case, the need arose for 
neutron source with neutron output the order 107 n/sec, which 
would require about 100 plutonium. Be/Pu atom ratio about 18:1 
was estimated for the completed source, which had neutron output 
n/sec corresponding 37.4 n/10° alphas 92% the theoretical neutron 
yield. These results confirmed the applicability the small-scale work the 
preparation larger sources. detailed account the fabrication this 
source has been published Wauchope (1959). 


CONCLUSIONS 

Neutron sources can prepared reducing with and sintering 
for minutes 1450 1250° depending the concentration. 

suitable size particle for preparation sources with high neutron 
yield was —200 mesh. 

Pressing the reagents 20,000 p.s.i. before sintering simplified the method 
eliminating separate degassing operation. 

The sources, generally prepared under vacuum, have also been prepared 
tank argon with neutron yields lower than equivalent sources sintered 
vacuum. 

The optimum neutron yield per unit volume sinter produced 1450° 
was obtained Be/Pu atom ratio the region 20:1. 


_ 
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NOTES 


THE SOURCE NOISE EBULLIOMETRY 


BLACKMORE 


INTRODUCTION 


part study ebulliometers (Blackmore 1959), was necessary 
determine whether the noise observed was produced the ebulliometers 
themselves the temperature-sensing elements. Since these were two 
thermistors Wheatstone bridge circuit there was some suspicion that the 
observed noise might thermistor ‘‘current (Brophy 1954). deter- 
mine which source was responsible, first dummy wire-wound bridge was 
tested different amplifier gain settings and different bridge currents. Then 
the thermistor bridge was tested the same bridge currents and amplifier gain 
settings with the thermistors instantaneously maintained the same temper- 
ature special furnace. Finally the noise was measured with the thermistor 
bridge the ebulliometers themselves and found significantly greater 
there. 


APPARATUS 

d-c. Wheatstone bridge (Tinsley Co. Type 3352) was used conjunction 
with copper-block furnace (Blackmore, Part 1959) calibrate the ther- 
mistors measuring their resistances different temperatures. The 
Wheatstone bridge unbalance voltages were detected with 
(Model 9835-B) d-c. microvolt amplifier and displayed Brown 
recorder (10 full scale). The temperature the copper-block furnace was 
measured with platinum resistance thermometer and Mueller bridge. 

After the thermistors V.E. Co.) had been calibrated matched pair 
was used a-c. bridge. The other two arms the bridge were made 
manganin wire accurately known resistance, about equal the resistance 
the thermistors the operating temperature, and were mounted 
thermally and electrically insulated can. The balancing decade boxes were 
Tinsley Co. (Type C-126) with 10X10, and 0.01 ohm steps. The boxes 
were oil-filled for temperature stability. Tinsley Co. (Type C-133 decade 
box with steps and 0.0001 ohms was also used. All leads were 
shielded and all shields were connected one ground point avoid ground 
loops. 

The circuit was arranged that either the thermistor bridge wire- 
wound dummy bridge could tested. The dummy bridge consisted four 
manganin wire arms wound heavy brass former screwed into heavy brass 
canister. This insured that all four arms were the same temperature. The 
equivalent resistance the bridge was about 150 ohms, which was the value 
expected for the thermistor resistances the operating temperature. 
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The a-c. waveform for the bridge was produced double-pole double- 
throw chopper series with 2-volt battery and current-limiting 
resistor. The choppers were manufactured the Tinsley Co. the design 
Dauphinee (1953). The bridge output was fed into Liston—Becker Model 
breaker amplifier (with its own choppers disconnected). This was done because 
other a-c. amplifier was conveniently available the laboratory that 
time. The amplifier output was then synchronously chopped and fed through 
filter network Brown Electronik recorder. The 
choppers were operated cycles because the amplifier designed operate 
this frequency. 

RESULTS 

The dummy-bridge sensitivity (Table chart divisions per ohm change 

the balancing arm shown function bridge current for different 


amplifier gain settings. The figures beside each entry the table are measure- 
ments the peak-to-peak noise that setting. 


TABLE 
Dummy-bridge sensitivity (div/Q) 


Bridge current 


gain setting 100 250 500 
16-1 2,950+0 14,800+1/4 
20-1 24,500+1 112,000+2 


test the thermistor bridge both thermistors were frozen into solder 
plug (50:50 m.p. 225° side side along with platinum resistance 
thermometer. This plug was then submerged oil the copper-block furnace 
and maintained temperatures comparable the operating temperatures 
expected the usual ebulliometer. Once again the noise was measured 
different amplifier gain settings and different bridge currents. will 
observed (Table that this noise was slightly larger than had been with 
the dummy bridge. This extra noise probably current 
(Brophy 1954). were extra pickup then the noise lower currents and 
amplifier gains should have been greater for the thermistor bridge than the 
dummy bridge. because their position the furnace the thermistors 
were more exposed hum pickup then turning the heater power off the 
furnace would have changed the noise the trace. This did not happen. 


TABLE 


Thermistor bridge sensitivity (div/Q) 


Bridge current 


Amplitier 
gain setting 100 250 500 
14-1 3,100+1/4 6,200 
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convert the noise measurement for the thermistor bridge equivalent 
temperature fluctuation quite simple. will treat the noise +1/2 
division gain 16-1 and bridge current though were all produced 
temperature fluctuation only one thermistor the bridge. 8000 
divisions per ohm (from Table change +1/2 division represents 
peak-to-peak noise +62.5 (micro-ohms). Now from the thermistor 
calibration known that this has temperature coefficient about 
100° Therefore the noise +62.5 equivalent +15 peak- 
to-peak. But know (Blackmore 1959, ref. that even for pure solvent 
typical ebulliometer the noise level about +150 
conclude therefore that the noise truly produced the ebulliometer and 
not the electrical equipment the thermistors. 
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SOME SENSITIVITIES ION GAUGES 
McGowan AND LARKIN KERWIN 


Schulz (1957) has published some characteristics the Bayard—Alpert 
ionization gauge used pressures above Hg. this pressure region, 
has been the practice use ionization gauges such the RCA 1949. has 
become known that many triodes are suitable for such measurements. 
particular, the 826-A twin-filament triode which currently available the 
surplus market about one dollar has been used several laboratories. 


PROCEDURE 

Two 826-A triodes, type gauge (Veeco RG-75), and 
standard Macleod gauge were exposed the same atmospheres various 
gases. recommended Schulz, the tubes were operated 100 electron 
current instead the usual order permit the application high 
pressure. 

With the plates (ion collectors) the tubes kept ground potential, the 
filaments were biased +22.5 volts and the grids +150 volts. Pressures 
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were varied from Hg. Ion currents from amp 
were measured with G-R model 1230-A electrometer. 

The following considerations are discussed length Schulz (1957). For 
given gas and given electron current, the number ions formed varies 
the pressure. The variation essentially linear one about 0.1 Hg. 
This has been verified the case the 826-A’s 0.01 Hg. For given 
pressure, the ion current varies linearly with the electron current satura- 
tion. For given pressure and electron current, different gases will give rise 
different ion currents. The ratio such ion currents for two gases should 
the same the ratio their cross sections for ion formation when bom- 
barded electrons known energy. shall refer the latter 
ence ratio” although measured one. Such cross sections are given 
Massey and Burhop (1952). The electron energy involved varies our case 
from 127.5 about the ionization potential, depending where between 
the grid and collector the ions are formed. For some gases the ratio cross 
sections moderately constant over much this range. 

From these considerations, the ion current seen proportional the 
electron current and the pressure. The sensitivity the gauge may then 
defined from the relation 


14 1_pS 


where may given units (mm For comparison purposes, argon 
may used standard, when the ratio (of the sensitivity for the gas 
that for argon) may compared with the reference ratio where 
the cross section for ion formation. 


RESULTS 
TABLE 


Observed and calculated sensitivities and ratios 


Sensitivity (mm Ratio 


Ref. cross This note This note 

Gas 1952) (100 ev) B-A* ratio 
Helium 0.35 2.5 0.5s 0.35 0.10 0.14 0.16 0.13 
2.87 12.0 2.8 0.82 0.67 0.80 0.75 
2.89 2.5 0.82 1.0 0.94 
Air 2.87 3.3 0.82 0.94 


*B-A: Bayard-Alpert. 


this work, have used the value cm? for the apparent cross 
section for ion formation argon. Schulz appears have used the value 2.64 
corresponding singly charged ions. The value 3.5 brings all 
published ratios into better agreement with the reference ratios except for He. 
high value for helium also evident our work. There little doubt that 
this due the presence impurities, very small amount which would 
have considerable effect the case helium. The other measured ratios are 


= 
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good agreement with the reference ratios when the experimental errors both 
this work and that used obtain the values the cross sections are 
taken into account. 

The sensitivity the Veeco gauge published here much lower than the 
official value (10 for fact, was found that two series sensitivities 
could found for this gauge when used high pressures. The sets values 
were consistent with cross sections, and the tube could from 
high low ion current positioning magnet near the tube, and other 
means which appeared suggest two stable configurations space 
charge the tube. This effect has been noted elsewhere (Carmichael, private 
Carter 1959, and references cited therein) and has been claimed 
due rather the charging the glass envelope either filament 
plate potential. have reported here the low values which were the more 
stable the range pressures used us. The 826-A tube did not display 
this dual sensitivity. 


Gratefully acknowledged assistance the measurements was given Paul 
Marmet and Jean-Denis Carette. The work part program carried out 
with the assistance N.R.C. Grant No. BT-45. 
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THE STATUS THE CANADIAN NEUTRON STANDARD* 
GEIGER 


The neutron emission rate the radium beryllium source No. N-200-1 
the National Research Council has recently been measured absolutely 
two methods (Geiger and Whyte 1959; Geiger 1959). The source was pre- 
pared 1948 and contained 183 radium element the form the 
carbonate, mixed with 2000 beryllium powder and pressed into pellet 
ensure physical stability. The pellet contained cylindrical Monel 
capsule 2-mm wall thickness, both the diameter and the height being 
mm. 

The first method calibration, based slowing down and capturing the 
neutrons July 1958. For the second method auxiliary polonium 
fluorine source was used and following the growth the radioactive 

*Issued No. 5598. 
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reaction product, sodium-22, the neutron emission rate could established. 
This source was then compared, manganese sulphate bath, with the radium 
beryllium source with resulting emission rate (3.16+0.10) 
the same date. Since the spectra radium beryllium source and 
polonium fluorine source are very different (Szilvasi al. 1960) 
somewhat larger error had assigned this value because uncertainties 
the intercomparison. The weighted mean both results is: 


Calibrations carried out different laboratories might expected give 
somewhat different values for the neutron emission rate and comparisons 
between various laboratories are therefore desirable. This facilitated the 
fact that relative intercomparisons sources can carried out with greater 
accuracy than absolute calibration. Therefore, only necessary com- 
pare the absolute standards the laboratories concerned. Probably the first 
results such intercomparisons (in 1951 and 1954) have been published 
Hughes (1954). Whereas the maximum spread calibrations was 20% 
1951 was reduced 10% 1954, showing that definite improvements 
neutron source calibrations were made. The most recent set intercomparisons 
was carried out Larsson (1958), who had compared the Swedish standard 
with standards from U.S.S.R., Belgium, Switzerland, U.S.A., and United 
Kingdom. The maximum spread between calibrations found was 4%. 

During the first half 1959 the Canadian neutron standard was compared 
with standards various laboratories. The prime requirement when such 
measurements are carried out make sure that the emission rate the 
Canadian standard did not change during the course these intercomparisons, 
possibly caused rough handling shipment. Measurements carried out 
this laboratory after completion the intercomparisons indicated that 
change neutron emission rate occurred within the precision the measure- 
ment +0.3%. 


RELATIVE SOURCE STRENGTH 


COUNTRY AND YEAR 
ABSOLUTE CALIBRATION 


SWEDEN|BELGIUMIGERMAN U.K. U.S.A. [CANADA 
1953-54] 1952 1958 1959 1955 1958 

Fic. Relative source strengths neutron standards various laboratories relation 
the Canadian standard (arbitrarily set equal 1.000). 
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The results the intercomparisons are shown Table and Fig. Column 
this table gives the laboratory where standard was calibrated absolutely 
and with which the Canadian standard was compared. 


TABLE 
Absolute source 
Date and place strength at date of Ratio of 
of inter- intercomparison (x)scale to 
Absolute neutron comparison Measured —_— ————————— QN RC. N.R.C. scale 
standard and year with N.R.C. ratio, Qz, ON.RC., on (x)scale, (relative source 
(1) (2) (3) (4) (5) (6) (7) 
Sweden Stockholm 0.843 2.687 3.214 3.187 0.992 
AB Atomenergi Jan. 1959 +0.8% +2% 
1953-1954 
Belgium Brussels 2.485 8.057 3.215 3.242 1.008 
Union Miniére du Feb. 1959 +0.12% +2% 
Haut Katanga 
1952 
Germany Braunschweig 0.6038 1.954 3.216 3.236 1.006 
Physikalisch Technische Mar. 1959 +0.3% +2% 
Bundesanstalt 
1958 
U. K. Teddington 3.215 3.128 0.972 
National Physical June 1959 +1.5% 
Laboratory 
U. S.A. Ottawa 0.3654 1.187 3.220 3.248 1.009 
National Bureau of July 1959 +1% +2% 
Standards 
1955 
Canada 3.208 3.205 1.000 
National Research (July) 
Council 
1958 


Column shows where the intercomparisons were carried out and column 
shows the ratio found. All except the U.S. source were radium beryllium 
sources which have similar neutron spectra. This allows intercomparisons 
made without great difficulties. The U.S. source radium beryllium 
source with consequently much less energetic neutron spectrum. The 
comparison was out using manganese sulphate bath activation and the 
same corrections for fast neutron losses from the source were applied 
described previous paper (Geiger and Whyte 1959). 

Column shows the values obtained for the absolute calibrations, corrected 
for the date intercomparison. The Swedish source was calibrated Larsson 
(1954, 1955) using four different methods: associated particle counting from 
the reaction, absolute counting with proton recoil chamber, 
and neutron density measurements aqueous solutions with small BF; 
chambers and gold foils. The Belgian source was calibrated Troyer and 
Tavernier (1954) neutron density measurement water with gold foils. 
The emission rate the German source was determined Kolb and von 
Droste (1959) using the same method. Since British standard was not 


— 
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available the time, Axton (1959) calibrated the Canadian source directly 
activation manganese sulphate bath. The U.S. source N.B.S. was also 
calibrated the bath method and, addition, neutron density measure- 
ment water using foils (De Juren al. 1955a, b). N.B.S. was very carefully 
compared with N.B.S. (Mosburg 1959). 

Column gives the neutron emission rate the Canadian standard corrected 
for the increase neutron emission due growth polonium-210 within the 
source and column shows its emission rate found comparison with the 
other standards. 

The last column shows the relative neutron emission rates for the Canadian 
standard (relative source strength) obtained the various laboratories; 
the Canadian value has arbitrarily been set equal 1.000. Figure shows this 
graphical form with only the errors the absolute calibration given. 

The agreement between the calibrations satisfactory, the maximum 
spread being 3.7%. Although, considering this relatively small spread, the 
errors quoted the various authors appear rather large has 
kept mind that calibrating neutron sources certain corrections are 
applied which are identical some the calibration methods. error this 
correction will therefore produce systematic shift all such calibrations. 
This applies particular the correction for fast neutron capture oxygen 
the case bath methods. Any further improvements the standardization 
neutron sources depends mainly the better knowledge such correction 
factors. 
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HIGH-ENERGY GAMMA RADIATION FROM 


Following recent coincidence study 253-day (Taylor and Scott 
1959) search was made for weak high-energy gamma-ray component 
observed der Mateosian and Goldhaber (1950), Wilson (1950), and 
Dzhelepov and (1956) using photodisintegration techniques. The 
energy the radiation determined these workers was found between 
1.66 and 2.2 Mev. intensity photons per disintegration was 
quoted Wilson, and photons per disintegration Dzhelepov and 

The observation gamma radiation from with scintillation counter 
this energy region complicated the presence prominent sum lines 
the gamma-ray pulse height spectrum 1596, 1825, 2041, 2166, and 2270 kev 
(Taylor and Scott 1959). Since the single counting rate proportional the 
solid angle and the sum coincidence rate the present investigation was 
carried out with relatively small solid angle. The detector consisted 
in. diam. in. long Nal crystal mounted Mont 6292 photo- 
multiplier tube. was mounted vertically lead castle with walls in. 
thick and open top. The source, 0.4 formed slow-neutron 
irradiation spec. pure silver, was placed 175 vertically above the crystal. 
The pulses from the counter were amplified with linear amplifier (resolving 
time and analyzed with Baird-Atomic 20-channel pulse height 
analyzer. 

The pulse height distribution (corrected for background pulses) the range 
1300 1800 kev shown Fig. The shapes the photopeaks (dashed 
lines) 1510 and 1560 kev were calculated using the gamma-ray intensities 
Dzhelepov (1957, 1958). The shapes the sum lines 1596 and 1825 kev 
were calculated from geometrical considerations and the known intensities 
the contributing low-energy gamma-rays (656, 885, and 940 kev). Compton 
distributions have been excluded for the sake simplicity since they not 
affect the shape the spectrum substantially 1700 kev. The solid curve 
the figure represents the sum the photopeaks and sum lines mentioned above. 
good agreement with the experimental points only one includes 
additional peak kev, shown. (Of course, the spectrum obtained 
using photopeaks and sum lines only will not exhibit the smearing due 
Compton distributions; this evident from the figure.) The gamma ray 
associated with this feature, assuming not spurious, has intensity 
(6+2) photons per disintegration, value which agrees with the observa- 
tions Wilson, and Dzhelepov and 

About one half-life after the first experimental runs, this region the 
spectrum was re-examined with the same geometry, etc. The counting rates 
the maxima the 1385- and 1510-kev photopeaks were reduced 
their original values, expected. The corresponding reduction for the 
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COUNTS PER MINUTE 


1300 1400 1600 1700 
GAMMA-RAY ENERGY (KEV) 


Pulse height spectrum between 1300 and 1800 obtained with 
Source-to-crystal distance, 175 em. 


peak was the error being due solely the low 
counting rate. seems improbable that the peak 1740 kev due 
accidental summing pulses the amplifier since this would result 
reduction factor The only accidental sum peaks expected this region 
the spectrum have energies 1700 and 1755 kev. 
Their intensities have been estimated less than 10% the 
peak. the other hand, the only genuine sum this region are those 
and J825 with intensities shown Thus would appear 
that the peak 1730 not due either genuine accidental super- 
position pulses due the strong gamma rays present. 
seems unlikely due contaminant view the high purity the 


silver used. 
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The peak 1730 kev does not correspond the crossover transition 
energy 1.8 Mev mentioned Nuclear Data Cards. Furthermore, this feature 
cannot the Compton maximum for gamma ray this energy. Taking the 
inflection point 1770 kev gives gamma-ray energy about Mev. 
check this possibility, however, the region above 1800 kev was examined 
with the spectrometer. Although the counting rate was found slightly 
excess the background, structure was observed the pulse height 
spectrum which had intensity appropriate Compton maximum. 
The counting rates encountered were comparable that 1800 kev and are 
probably due weak genuine sum lines with energies excess 1825 kev 
(Taylor and Scott 1959). 

The observations reported here suggest that the observed peak due 
single gamma ray energy 1730+30 kev emitted Since its energy 
lies the range 1.6 2.2 Mev, and gamma ray comparable intensity 
observed 1825 kev, probably the same radiation reported other 
workers. 

transition this energy cannot included the level scheme 
the present time. interest note, however, that other weak gamma rays 
energies 575, 541, 499, 471, and 437 kev have been reported Cork 
(1950), none which have been incorporated into the level scheme. Two 
these gamma-ray energies, when added 1730 kev, give sum energies worthy 
note: 1730+541 2271 kev and 1730+437 2167 kev. Sum lines these 
energies previously reported the authors are known arise from genuine 
coincidences between intense gamma rays: 656+1510 2166 kev, 
760+1510 2270 kev. the 1730- and gamma rays 
cascade, there must level 1730 kev (see Taylor and Scott (1959) 
for level scheme). Similarly, cascade containing the and 
gamma radiations implies level kev. (The alternative, level 
1197 kev, would conflict with the systematic trends the energy levels 
nuclei.) The inclusion other weak gamma rays into the decay 
scheme not facilitated either possibility, however, and the clarification 
the situation must await further work with the weak gamma radiations. 
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Frequency Measurement Standard Frequency 


Measurements are made Ottawa, Canada, using N.R.C. caesium-beam frequency resona- 
tor reference standard (with assumed frequency 192 631 770 c.p.s.). Frequency 
deviations from nominal are quoted parts per negative sign indicates that the fre- 
quency below nominal. 


GBR, 16 kc/s 
Date, 
January MSF, 84-hour 24-hour KK2XEI, 
1960 60 kce/s average* average 60 kc/s 
1 — 167 -171 —169 N.M. 
2 -—173 —170 —170 N.M. 
3 -—170 -171 — 169 N.M. 
4 —173 —168 —169 —83 
5 — 169 —168 — 167 —82 
6 -173 -170 — 84 
7 — 167 —172 -171 —85 
8 —175 —168 -—170 —85 
9 -—173 —170 —172 N.M. 
10 —173 —170 —170 N.M. 
11 —167 —169 —169 — 84 
12 —165 -171 -171 —87 
13 —172 —170 — 169 — 85 
14 —166 —170 -—172 —85 
15 —176 —172 —170 —85 
16 —170 —170 N.M. 
17 — 168 —173 —172 N.M. 
18 —188 -171 —170 —87 
19 N.M. -—172 —170 —83 
20 N.M. —167 — 168 —84 
21 N.M. —-175 —172 —88 
22 N.M. —170 -171 N.M. 
23 —176 —174 —175 N.M. 
24 -177 —178 —176 N.M. 
25 -173 —178 —176 N.M. 
26 —175 —173 —174 N.M. 
27 —176 —175 —175 N.M. 
28 —180 -177 —174 —89 
29 -—179 —174 —174 —89 
30 —175 N.M. 
31 — 182 -177 -177 N.M. 
Midmonthly mean —174 —172 —172 —85 


Midmonthly mean 
of WWV — 159 


Note: N.M. no measurement. 
*Time of observations: 00.00 to 07.00 and 22.30 to 24.00 U.T. 
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